Search for: All records

Context.One of the most important open questions in planet formation is how dust grains in a protoplanetary disk manage to overcome growth barriers and form the ∼100 km planet building blocks that we call planetesimals. There appears to be a gap between the largest grains that can be produce by coagulation, and the smallest grains that are needed for the streaming instability (SI) to form planetesimals. Aims.Here we explore the novel hypothesis that dust coagulation and the SI work in tandem; in other words, they form a feedback loop where each one boosts the action of the other to bridge the gap between dust grains and planetesimals. Methods.We developed a semi-analytical model of dust concentration due to the SI, and an analytic model of how the SI affects the fragmentation and radial drift barriers. We then combined them to model our proposed feedback loop. Results.In the fragmentation-limited regime, we find a powerful synergy between the SI and dust growth that drastically increases both grain sizes and densities. We find that a midplane dust-to-gas ratio ofϵ ≥ 0.3 is a sufficient condition for the feedback loop to reach the planetesimal-forming region for turbulence values 10⁻⁴ ≤ α ≤ 10⁻³and grain sizes 0.01 ≤ St ≤ 0.1. In contrast, the drift-limited regime only shows grain growth without significant dust accumulation. In other words, planetesimal formation remains challenging in the drift-dominated regime and dust traps may be required to allow planet formation at wide orbital distances. 

Abstract The streaming instability (SI) is a leading mechanism for concentrating solid particles into regions dense enough to form planetesimals. Its efficiency in clumping particles depends primarily on the dimensionless stopping time (τ_s, a proxy for particle size) and dust-to-gas surface density ratio (Z). Previous simulations identified a criticalZ(Z_crit) above which strong clumping occurs, where particle densities exceed the Hill density (thus satisfying a condition for gravitational collapse), over a wide range ofτ_s. These works found that, forτ_s≤ 0.01,Z_critwas above the interstellar medium value (∼0.01). In this work, we reexamine the clumping threshold using 2D axisymmetric, stratified simulations at high resolution and with relatively large (compared to many previous simulations) domain sizes. Our main results are as follows: First, whenτ_s = 0.01, strong clumping occurs even atZ ≲ 0.01, lower thanZ_critfound in all previous studies. Consequently, we revise a previously published fit to theZ_critcurve to account for this updatedZ_crit. Second, higher resolution results in a thicker dust layer, which may result from other instabilities manifesting, such as the vertically shearing SI. Third, despite this thicker layer, higher resolution can lead to strong clumping even with a lower midplane dust-to-gas density ratios (which results from the thicker particle layer) so long asZ ≳ Z_crit. Our results demonstrate the efficiency of the SI in clumping small particles atZ ∼ 0.01, which is a significant refinement of the conditions for planetesimal formation by the SI. 

Abstract In the theory of protoplanetary disk turbulence, a widely adopted ansatz, or assumption, is that the turnover frequency of the largest turbulent eddy, Ω_L, is the local Keplerian frequency Ω_K. In terms of the standard dimensionless Shakura–Sunyaevαparameter that quantifies turbulent viscosity or diffusivity, this assumption leads to characteristic length and velocity scales given respectively by $\sqrt{\alpha }H$and $\sqrt{\alpha }c$, in whichHandcare the local gas scale height and sound speed. However, this assumption is not applicable in cases when turbulence is forced numerically or driven by some natural processes such as vertical shear instability. Here, we explore the more general case where Ω_L≥ Ω_Kand show that, under these conditions, the characteristic length and velocity scales are respectively $\sqrt{\alpha /R\prime }H$and $\sqrt{\alpha R\prime }c$, where $R\prime \equiv {\mathrm{\Omega }}_L/{\mathrm{\Omega }}_K$is twice the Rossby number. It follows that $\alpha =\tilde{\alpha }/R\prime$, where $\sqrt{\tilde{\alpha }}c$is the root-mean-square average of the turbulent velocities. Properly allowing for this effect naturally explains the reduced particle scale heights produced in shearing box simulations of particles in forced turbulence, and it may help with interpreting recent edge-on disk observations; more general implications for observations are also presented. For $R\prime >1$, the effective particle Stokes numbers are increased, which has implications for particle collision dynamics and growth, as well as for planetesimal formation. 

Abstract Disk vortices, seen in numerical simulations of protoplanetary disks and found observationally in Atacama Large Millimeter/submillimeter Array and Very Large Array images of these objects, are promising sites for planet formation given their pebble trapping abilities. Previous works have shown a strong concentration of pebbles in vortices, but gravitational collapse has only been shown in low-resolution, two-dimensional, global models. In this Letter, we aim to study the pebble concentration and gravitational collapse of pebble clouds in vortices via high-resolution, three-dimensional, local models. We performed simulations of the dynamics of gas and solids in a local shearing box where the gas is subject to convective overstability, generating a persistent giant vortex. We find that the vortex produces objects of Moon and Mars mass, with a mass function of power-law $d\ln N/d\ln M=-1.6\pm 0.3$. The protoplanets grow rapidly, doubling in mass in about five orbits, following pebble accretion rates. The mass range and mass doubling rate are in broad agreement with previous low-resolution global models. We conclude that Mars-mass planetary embryos are the natural outcome of planet formation inside the disk vortices seen in millimeter and radio images of protoplanetary disks. 

Abstract Given the important role turbulence plays in the settling and growth of dust grains in protoplanetary disks, it is crucial that we determine whether these disks are turbulent and to what extent. Protoplanetary disks are weakly ionized near the midplane, which has led to a paradigm in which largely laminar magnetic field structures prevail deeper in the disk, with angular momentum being transported via magnetically launched winds. Yet, there has been little exploration of the precise behavior of the gas within the bulk of the disk. We carry out 3D, local shearing box simulations that include all three low-ionization effects (ohmic diffusion, ambipolar diffusion, and the Hall effect) to probe the nature of magnetically driven gas dynamics 1–30 au from the central star. We find that gas turbulence can persist with a generous yet physically motivated ionization prescription (order unity Elsässer numbers). The gas velocity fluctuations range from 0.03 to 0.09 of the sound speedc_sat the disk midplane to ∼c_snear the disk surface, and are dependent on the initial magnetic field strength. However, the turbulent velocities do not appear to be strongly dependent on the field polarity, and thus appear to be insensitive to the Hall effect. The midplane turbulence has the potential to drive dust grains to collision velocities exceeding their fragmentation limit, and likely reduces the efficacy of particle clumping in the midplane, though it remains to be seen if this level of turbulence persists in disks with lower ionization levels. 

Abstract The streaming instability (SI) is a leading candidate for planetesimal formation, which can concentrate solids through two-way aerodynamic interactions with the gas. The resulting concentrations can become sufficiently dense to collapse under particle self-gravity, forming planetesimals. Previous studies have carried out large parameter surveys to establish the critical particle to gas surface density ratio (Z), above which SI-induced concentration triggers planetesimal formation. The thresholdZdepends on the dimensionless stopping time (τ_s, a proxy for dust size). However, these studies neglected both particle self-gravity and external turbulence. Here, we perform 3D stratified shearing box simulations with both particle self-gravity and turbulent forcing, which we characterize via a turbulent diffusion parameter,α_D. We find that forced turbulence, at amplitudes plausibly present in some protoplanetary disks, can increase the thresholdZby up to an order of magnitude. For example, forτ_s= 0.01, planetesimal formation occurs whenZ≳ 0.06, ≳0.1, and ≳0.2 atα_D= 10⁻⁴, 10^−3.5, and 10⁻³, respectively. We provide a single fit to the criticalZrequired for the SI to work as a function ofα_Dandτ_s(although limited to the rangeτ_s= 0.01–0.1). Our simulations also show that planetesimal formation requires a mid-plane particle-to-gas density ratio that exceeds unity, with the critical value being largely insensitive toα_D. Finally, we provide an estimation of particle scale height that accounts for both particle feedback and external turbulence. 

Abstract Kuiper Belt objects (KBOs) show an unexpected trend, whereby large bodies have increasingly higher densities, up to five times greater than their smaller counterparts. Current explanations for this trend assume formation at constant composition, with the increasing density resulting from gravitational compaction. However, this scenario poses a timing problem to avoid early melting by decay of²⁶Al. We aim to explain the density trend in the context of streaming instability and pebble accretion. Small pebbles experience lofting into the atmosphere of the disk, being exposed to UV and partially losing their ice via desorption. Conversely, larger pebbles are shielded and remain icier. We use a shearing box model including gas and solids, the latter split into ices and silicate pebbles. Self-gravity is included, allowing dense clumps to collapse into planetesimals. We find that the streaming instability leads to the formation of mostly icy planetesimals, albeit with an unexpected trend that the lighter ones are more silicate-rich than the heavier ones. We feed the resulting planetesimals into a pebble accretion integrator with a continuous size distribution, finding that they undergo drastic changes in composition as they preferentially accrete silicate pebbles. The density and masses of large KBOs are best reproduced if they form between 15 and 22 au. Our solution avoids the timing problem because the first planetesimals are primarily icy and²⁶Al is mostly incorporated in the slow phase of silicate pebble accretion. Our results lend further credibility to the streaming instability and pebble accretion as formation and growth mechanisms. 

Abstract Pebble accretion is recognized as a significant accelerator of planet formation. Yet only formulae for single-sized (monodisperse) distribution have been derived in the literature. These can lead to significant underestimates for Bondi accretion, for which the best accreted pebble size may not be the one that dominates the mass distribution. We derive in this paper the polydisperse theory of pebble accretion. We consider a power-law distribution in pebble radius, and we find the resulting surface and volume number density distribution functions. We derive also the exact monodisperse analytical pebble accretion rate for which 3D accretion and 2D accretion are limits. In addition, we find analytical solutions to the polydisperse 2D Hill and 3D Bondi limits. We integrate the polydisperse pebble accretion numerically for the MRN distribution, finding a slight decrease (by an exact factor 3/7) in the Hill regime compared to the monodisperse case. In contrast, in the Bondi regime, we find accretion rates 1–2 orders of magnitude higher compared to monodisperse, also extending the onset of pebble accretion to 1–2 orders of magnitude lower in mass. We find megayear timescales, within the disk lifetime, for Bondi accretion on top of planetary seeds of masses 10⁻⁶to 10⁻⁴M_⊕, over a significant range of the parameter space. This mass range overlaps with the high-mass end of the planetesimal initial mass function, and thus pebble accretion is possible directly following formation by streaming instability. This alleviates the need for mutual planetesimal collisions as a major contribution to planetary growth. 

ABSTRACT Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination i < 90°) decrease inclination, as well as semimajor axis (a) and eccentricity (e) until orbital alignment with the gas disc (‘disc capture’). Captured stellar-origin black holes (sBH) add to the embedded AGN population that drives sBH–sBH mergers detectable in gravitational waves using LIGO–Virgo–KAGRA or sBH–SMBH mergers detectable with Laser Interferometer Space Antenna. Captured stars can be tidally disrupted by sBH or the SMBH or rapidly grow into massive ‘immortal’ stars. Here, we investigate the behaviour of polar and retrograde orbiters (i ≥ 90°) interacting with the disc. We show that retrograde stars are captured faster than prograde stars, flip to prograde orientation (i < 90°) during capture, and decrease a dramatically towards the SMBH. For sBH, we find a critical angle iret ∼ 113°, below which retrograde sBH decay towards embedded prograde orbits (i → 0°), while for io > iret sBH decay towards embedded retrograde orbits (i → 180°). sBH near polar orbits (i ∼ 90°) and stars on nearly embedded retrograde orbits (i ∼ 180°) show the greatest decreases in a. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial a. For sBH, disc capture time is longest for polar orbits, low-mass sBH, and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the spin distribution of embedded sBH.